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METHOD AND ARRANGEMENT FOR DETE RMINING AN OBJECT CONTOUR 

Cross reference to related application: This application claims benefit from Swedish patent 
application no SE-03021 14-4, filed July 21, 2003. 

Field of the Invention 

The present invention relates to a method of determining a sought object contour in a digital 
microscope image according to the preamble of claim 1. The invention also concerns a 
corresponding arrangement according to the preamble of claim 5 and a digital storage 
medium comprising a corresponding computer program according to the preamble of claim 6. 

Background Art 

The blood of himians contains red blood cells, platelets and white blood cells. In medical 
applications there is often an interest in classifying the white blood cells into approximately 
twenty different sub classes and in estimating the distribution of these subclasses for 
individual patients. Traditionally this classification may be performed in a laborious way by a 
human expert looking through the eyepiece of a brightfield microscope and thereby studying 
blood which has been wedged on a slide. The classification may also, with less separable sub 
classes though, be performed in a so called cell coimter where diluted blood passes, cell by 
cell, through a narrow strait in a so called flow cell. While the cells pass through the flow cell 
parameters like electrical impedance and light scattering are measured, which subsequently 
are the basis for the classification of the actual white blood cell performed by the cell counter. 

The advances within image analysis during recent years have made it possible to produce 
automatical brightfield microscopes that are capable of scaiming a blood smear for white 
blood cells, autofocus on them, store well focused cell images, perform automatical 
segmentation of the resulting cell images, calculate so called features for each segmented cell 
image and finally classify each found white blood cell based on the information in the 
calculated features. 

Segmentation is the process of splitting the contents of an image into different regions. In the 
case of white blood cells, the desired result is that white blood cells are separated from the 
other objects in the image like background, red blood cells, platelets, garbage and, of course. 
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from other white blood cells. It may also be desirable that the segmentation separates the 
nucleus of the white blood cell from its cytoplasm. 

The main reason for segmenting the individual cells and their parts is that there, today, do not 
exist any classifaction methods that work well enough and quickly enough without the 
separation of image information that the segmentation corresponds to. 

Acceptable segmentation results may sometimes be obtained by using so called binary 
thresholding, meaning that each image element, also called pixel, in the cell image is dealt 
with depending on if its intensity is greater than or below a threshold intensity. The threshold 
intensity may be local or global, i.e. constant over the whole image. 

An apparent drawback of binary thresholding is that the object contour that is obtained is 
sensitive to the level of the chosen threshold intensity, for example when the illumination 
conditions vary within av cell image. In addition, thresholding is performed regardless of 
shape or so called non-local information, i.e. how the object appears at some distance from 
the sought contour. Binary thresholding may therefore give the result that adjacent objects are 
not separated but segmented as one object or that an object with an original smooth, rounded 
contour gets a zigzag contour after segmentation. Binary thresholding may also give the result 
that a single original object is segmented as multiple objects. Binary thresholding and more 
sophisticated thresholding methods are therefore most often only used as initial operations in 
more powerful segmentation methods. 

So called active contour models is a class of methods that have bera used for quite a while for 
segmentation of inoiages. Using the active contour models the contour of the object is 
determined directly instead of indirectly via which image elements that belong to the object, 
like it is done using thresholding. All active contour models use some kind of contour that 
iteratively, during the process of segmentation, is altered under the influence of so called 
intemal forces that originate from the contour and its shape and under the influence of so 
called external forces that originate from the image and mainly from its contents of edge 
information. 

One way of enhancing edge information in images is to calculate so called gradient images. 
Briefly that means that elements in the image, where the intensity is increased or decreased 
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with a large amount on a distance corresponding to a few image elements, get 
correspondingly large values at the corresponding image element positions in the gradient 
image. Edges are a good example of objects that give rise to large gradient values. The 
gradient images may then be post-processed in order to get rid of isolated large values, that 
S probably were caused by garbage or small details, and in order to enhance long joint streaks 
of large values probably caused by edges of the sought object and other objects. An image 
where the values of the image elements consist of some kind of gradient values will be called 
an edge image below. 

10 The active models may be divided into two groups. The first group is the group of 
geometrical models like "Level Sets" and "Fast Marching" methods. The contour is 
represented by the set of image elements of the contour. The advantage of that representation 
is that the contour may have arbitrary shapes. The iterations comprise operations on the image 
elements of the contour and their corresponding closest neighbours. Since the contour may 

IS have arbitrary shapes, it is relatively complicated and computationally heavy to administrate 
all the image elements of the contour, their neighbours and their matual relative positions. 

The second group is the group of parametrical models like "snakes", see the paper "Snakes: 
Active contour models," Intemational Joumal of Computer Vision, 1(4): 321-331, 1987, by 

20 Kass, Witkin och Terzopoulos. A snake is a parametric contour model. The parameters in the 
model may for example consist of the positions of the comers of a polygon, but most often the 
model is more advanced, for example by consisting of joint curved contour elements. The 
internal forces of the snake strive for a smooth contour while the external forces strive for 
letting the contour pass close to msay of the image elements that contain edge information. 

25 During the process of segmentation, an algorithm based on snakes will itererate position and 
shape of the snake until a good compromise - hopefully the most fitting contour - between 
internal and external forces is found. By increasing the influence of the internal forces, the 
shape of the final contour estimate may become more or less limited to circular shaped 
objects. By correspondingly increasing the influence of the extemal forces, the shape of the 

30 final contom estimate may become more irregular. Since the snake is represented by a limited 
number of parameters, its possible shapes are limited. For white blood cells, which have fairly 
regular contours, this limitation is seldom a drawback. On the contrary, it is an advantage that 
the snake will have a regular shape even in those cases where the contour of the object is seen 
only vaguely in the cell image and where it therefore will be quite a distance between image 
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elements having edge infoimation. In addition, the limited number of parameters will lead to 
that the computational power needed during iterations will be less compared to that of a 
geometrical model. 

If the snake is started close to the sought contour, it will probably converge safely and with a 
small number of iterations. Such a start may be easy for a human operator but harder to 
accomplish using automatic image analysis. Automatic segmentation will therefore require 
that the snake may be started relatively far away from the sought object contour but still 
converge safely and quickly. That requirement will in turn lead to the requirement that the 
external forces of the snake must be able to lead the snake correctly even from image 
elements that are situated far away from the sought contour. In order to be useful as input to a 
snake, the external forces should simply, for each image element in the original image, be 
pointing towards image elements with a greater amount of edge information, i.e. with greater 
edge resemblance. Therefore it is suitable to use the expression "vector force field" for the 
external forces as a function of the location in the cell image. 

One previously known way of obtaining a vector force field is to start with an edge image 
with values f(x,y), that depend on the x and y coordinate of the corresponding image element 
and calculate differences in f(x,y) with respect to x and y. The calculations of differences may 
for exem|>le be defined by 

fjt(x,y)=f(x+l,yH(x-l,y) (eq. 1) and 

f_y(x,y)=f(x,y+l)-f(x,y-l) (eq. 2), respectively. 

The obtained vector force field [f jc(x,y),f_y(x,y)] is an example of an extemal force field. 
The greatest drawback of such a vector force field is that the amount of force rapidly descends 
as the coordinates (x,y) move away firom those points that have an edge image f(x,y) that 
differs from zero. Thus, such a vector force field does not work very well when the snake is 
started firom a distance relatively far away from the sought object contour. 

One previously known way of improving the ability of the vector force field to act at a 
distance is to smear the edge image, before or after the calulations of differences with for 
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example a two-dimensional Gaussian filter, but that may result in other problems like details 
that become less apparent. 

The Gradient Vector Flow Method is a previously know method that accomplishes a static 
extemal vector force field, a so called GVF field (Gradient Vector Flow Field), calculated on 
the basis of a, like shown above, difference calculated edge image contaming information on 
the sought object contour. (See for example the Master's thesis "Segmentation of 
Histopathological Tissue Sections Using Gradient Vector Flow Snakes", Centre for 
Mathematical Sciences, Lund University and Lund Institute of Technology, March 18 2002, 
by Adam Karlsson or the original reference "Gradient vector flow: A new external force for 
snakes," IEEE Proceedings on Computer Vision and Pattern Recognition, Puerto Rico, pages 
66-71, 1997 by Xu och Prince). 

A GVF-field fulfills the requirement of being an extemal vector force field that leads the 
snake in the correct direction also at large geometrical distances firom the sought contour. A 
snake that uses a GVF-field is called a GVF-snake below. The GVF-field is then calculated 
once per segmentation - at the beginning of it. Such segmentation methods are previously 
known. 

In the Master's thesis above by Adam Karlsson there is further described a fast way of 
iterating the parameters of the snake, i.e. how the position and the shape of the snake is 
iterated based on the GVF-field. In the thesis it is pointed out that the tinung performance of a 
GVF-snake, which uses this fast way of iterating the parameters, will be limited by the time it 
takes to calculate the GVF-field. 

The traditional way of calculating the GVF field according to Xu och Prince contains iterative 
solving of 

tiAu-(u-f.x)(f_x^+fj^)=0 (eq.3)and 



^lAv-(v-f.x)(fJ:2+f_y2)=0 (eq.4). 
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where A is the so called Laplace operator - a kind of two-dimensional second order difference 
operator - and where ^ is a parameter that may have to be tuned depending on the application, 
i.e* what kind of objects it is and in which environment they are situated. 

The calculation of the GVF-field means that a vector force field [f jc(x,y),f_y(x,y)] with 
limited ability to act at a distance results in another vector force field - the GVF-field, 
[u(x,y),v(x,y)] - with an improved ability to act at a distance without the loss of details. The 
input to the calculation can be regarded as two images, f^ and f jy respectively, and the 
output can be regarded as two other images, u and v respectively. 

Solving equation systems 3 and 4 requires a considerable amount of iterations in order to get 
a converged result If the area of interest in the original image has the size of m tunes n pixels, 
each of the equation systems 3 and 4, will contain m times n equations, which means just as 
many unknowns, which for images easily becomes ten thousand unknowns which shall be 
updated during each iteration. 

It is not unusual that the traditional way of calculating the GVF-field, with a IGHz PC- 
processor, will need several tenths of a second per cell image. For an analysis of 200 white 
blood cells that corresponds to a total segmentation time of close on one minute, which limits 
the performance of existing automatic brightfield microscopes unless extra, expensive 
hardware is added. 

Sunmiarv of the Invention 

An object of the present invention is to wholly or partly eliminate the above problentis. 

This object is achieved by a method of determining a sought object contour in a digital unage 
according to claim 1, by a corresponding arrangement according to claim 5 and a digital 
storage medium comprising a computer program according to claim 6. 

More specifically, the invention then relates, according to a first aspect, to a method of 
determining a sought object contour in a digital microscope image, which comprises a 
plurality of image elements and reproduces a biological material. The method is characterized 
by the steps of assigning edge values to at least a first subset of the image elements in the 
image; assigning values of a first gradient vector component whose values each comprise a 
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first linear combination of edge values of some surrounding image elements to at least a 
second subset of the image elements in the image; assigning values of a second gradient 
vector component whose values each comprise a second linear combination of edge values of 
some surrounding image elements to at least a third subset of the image elements in the 
image; and calculating an estimate of the sought object contour based upon values of the first 
and the second gradient vector components. 

Such a method has been found to result in a quickly computed approximative GVF field, 
denoted AGVF field below. The computation of the AGVF field is performed so that, 
although the resulting AGVF field differs somewhat from the traditional GVF field, the 
segmentation results for white blood cells appear to be mainly unaffected. Using a method 
according to the invention, the time for using snakes with GVF performance in the 
segmentation can thus be considerably lowered. The present invention therefore makes it 
possible to use an improved segmentation in commercially available automatical scanning 
microscopes for analysis of white blood cells. 

For other applications than white blood cells, AGVF may be used as it is or as a way of 
accelerating, for example by skipping the initial iterations in a traditional GVF method. 

According to a second aspect, the invention concerns an arrangement for determining a 
sought object contour in a digital microscope image, which comprises a plurality of image 
elements and reproduces a biological material. The arrangement is characterized by means for 
assigning edge values to at least a first subset of the image elements in the image; means for 
assigning values of a first gradient vector component whose values each comprise a first 
linear combination of edge values of some surrounding image elements to at least a second 
subset of the image elements in the image; means for assigning values of a second gradient 
vector component whose values each comprise a second linear combination of edge values of 
some surroimding image elements to at least a third subset of the image elements in the 
image; and means for calculating an estimate of the sought object contour based upon values 
of the first and the second gradient vector components. 

The arrangement gives advantages corresponding to those of the method and can also be 
varied similarly. 
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According to a third aspect^ the invention relates to a digital storage medium comprising a 
computer program for determining a sought object contour in a digital microscope image» 
which comprises a plurality of image elements and reproduces a biological material. The 
storage medivun is characterized by instructions corresponding to the steps of assigning edge 

5 values to at least a first subset of the image elements in the image; assigning values of a first 
gradient vector component whose values each comprise a first linear combination of edge 
values of some surrounding image elements to at least a second subset of the image elements 
in the image; assigning values of a second gradient vector component whose values each 
comprise a second linear combination of edge values of some surrounding image elements to 

10 at least a third subset of the image elements in the image; and calculating an estimate of the 
sought object contour based upon values of the first and the second gradient vector 
components. 

The computer program gives advantages corresponding to those of the method and can also 
IS be varied similarly. 

Brief Description of the Figures 

Figs la and lb show a digital gray scale image of a white blood cell surrounded by 
background, red blood cells and a platelet. The image is an example of what is meant by a cell 
20 image. 

Fig. 2 shows a preliminary thresholding of the nucleus of the white blood cell. 

Fig. 3 shows an edge image where the darkest areas are possible cell contours and the 
25 brightest areas are such that most likely are not cell contours. 

Fig. 4 shows, with arrows, one from the edge image obtained approximative GVF field 
superimposed on fig. 3 for comparison. 

30 Fig. 5 shows a possible seed contour for a snake - namely the so called convex hull of the 
preliminarily thresholded ceU nucleus of fig. 2. 

Fig. 6 shows the estimated object contour when the segmentation is finished with the 
cytoplasm, i.e. the cell contour. 
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Fig. 7 shows the estimated cell contour of fig. 6 together with the estimated contour of the cell 
nucleus. 

S Figs 8a and 8b show an arrangement according to the invention. 
Fig. 9 is a flow chart of a method according to the invention. 
Description of Preferred Embodiments 

10 

A white blood cell consists, from a segmentation point of view, of two parts - the cell nucleus 
22 and the surrounding cytoplasm 21. The results of segmentation of the two parts are to 
some amount dependent upon one another: 

15 In order to succeed m finding the border between cytoplasm and background using automatic 
image analysis in spite of the presence of adjacent cells, marked R and T in fig. lb, it is useful 
to be able to start the snake from a so called seed contour which is completely inside the cell. 
Therefore one wishes to have access to an estimate, for example a segmentation, of the cell 
nucleus as a so called seed contour for the snake. 

20 

In order to simplify the segmentation of the cell nucleus it is, on the other hand, good to have 
access to an inoage where there is only cytoplasm and cell nucleus left - i.e. an image where 
the cell already is segmented from the background and adjacent cells. To aviod an iterative 
process, a preliminary segmentation, see fig. 2, of the cell nucleus is used in one preferred 
25 embodiment. That segmentation is good enougih for its purpose. 

In one preferred embodiment a nimiber of partial operations are performed starting from an 
original cell image, which, like the image in fig. la, contains the white blood cell whose 
contour and nucleus are sought. Fig. 8b gives a survey of how means corresponding to these 
30 partial operations may relate to one another. 

PRELIMINARY SEGMENTATION OF THE NUCLEUS 

During the partial operation corresponding to 83, a preliminary segmentation of the nucleus is 
calculated in order to achieve a seed contour for a snake. Fig. 2 shows the result of such a 
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preliminary segmentation. This preliminary segmentation is performed using known binary 
thresholding methods preferrably operating on the green colour component of the original cell 
image or, possibly, on a weighted combination of green and blue colour component. 

CALCULATION OF THE EDGE IMAGE 

In partial operation corresponding to 84 the values of an edge image are calculated by for 
example a gradient operation followed by a processing that only lets the strongest edges 
remain in the image. See the darkest areas of fig. 3. 

In the edge image it is preferred to incorporate the preliminary nucleus as an "anti*'-edge with 
negative values of its edge image elements - see the brightest areas in fig. 3. By doing so, the 
edge image will, when the AGVF-field has been calculated, have added the information that 
the sought cell contour is situated outside the cell nucleus. 

In a preferred embodiment, the edge image is trinary, i.e. its elements have three possible 
values, one value for edge (1), one value for nucleus (-1) and one value for other (0). It is also 
possible to let the edge values and the nucleus values differ individually from the value for 
other (0). 

It is not necesssary to calculate the contents of the edge image for all image elements of the 
cell imtage. It is sufficient to calculate it in the surroundings of the preliminarily segmented 
cell nucleus, where the surroundings are large enough to safely contain a white blood cell. By 
not working with an uimecessarily large edge image, there is a cut in the amount of 
calculations needed - especially during the calculation of the AGVF field. 

It is preferrable to choose a quadratic edge image with side lenghts that are a potence of 2 
image elements, since the Fourier transforms below are most efficiently calculated for such 
side lenghts using the well known Fast Fourier Transform, FFT. Such a quadratic edge image 
is an example of a first subset of image elements. 

DETERMINATION OF SNAKE SEED 

In the partial operation corresponding to 85 a calculation of a seed contour is performed using 
data from a partial operation correspondning to 83 as the input. For segmentation of white 
blood cells the so called convex hull of the preliminarily segmented nucleus is preferrably 
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used. By the expression "convex hull** a set of image elements with a convex contour 
containing all the image elements of the preliminarily segmented cell nucleus is intended. See 
fig. 5. 

CALCULATION OF THE COMPONENTS OF THE AGVF FIELD 
In the partial operation 86 values for the components of the AGVF field are calculated using 
the edge image from a partial operation corresponding to 84 as input The calculations may be 
performed using at least two methods corresponding to one principle each. 

According to the first method, the edge hnage is differentiated first, which is preferrably done 
according to equations 1 and 2 above. 

Then the main AGVF calculation is performed. According to the invention, the AGVF field is 
calculated by solving 

fiAu-(u-f_x)-l=0 (eq. 5) and 

KiAv-(v.f_y) 1=0 (eq. 6), 

where A is the so called Laplace operator - a kind of two-dimensional second order difference 
operator, which preferrably is calculated by two-dimensional filtering with the 3 by 3 matrix 

0 1 0 
L = 1 -4 1 
0 1 0. 

Equations S and 6 may then be rewritten as 

(Hl^l)*u=-f ji (eq. 7) and 

(Ml-l)*v=-fjr (eq.8), 

where the symbol * means so called convolution of the 3 by 3 matrix (tiL-1) and the m by n 
images u and v respectively. The contents of the 3 by 3 matrix is thus a weighted combination 
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of a Laplace filter and a unity filter. This rewriting is the key to the speed in the present 
AGVF method. It is a well known fact that convolution is replaced by element wise 
multiplication during Fourier transform and vice versa. 

Such rewriting is not possible with the traditional GVF method due to that the factor (f_x^+ 
f-J^)» which varies over the image, exists in equations 3 and 4. 

Equations 7 and 8 may be quickly solved using the Fast Fourier Transform, FFT. The 
components u and v of the AGVF field may in that way, explicitly, be obtained through the 
expressions 

u = - F*[F(f^x)./F(nI^l)] (eq. 9) and 
V = . F*[F(f jr)./F(nI^l)] (eq. 10). 

where F denotes the two-dimensional discrete Fourier Transform, denotes the 
corresponding inverse Fourier Transform and ./ denotes element wise division in the 
frequency plane. 

As long as p. is constant, neither F(p.L-l) nor 1 ./ F(^L-l) will vary, which leads to that the 
whole partial result IJ F(pL-l) may be stored in order to speed up the calculations even 
more. 

Suitable values of \i may be in the range 10-30. Which values that are the best depends on the 
sought objects, the magnification of the objective, the cell size, the pixel size of the camera, 
etc. It is recommended to choose |x by experimenting with a nxmiber of cells under manual 
supervision of the segmentation. 

As can be seen in equations 9 and 10, the calculations of u and v are linear operations on f_x 
and f_y respectively, which in their turn come from linear differentiation operations on f . 
Therefore it is possible, according to a second possible calcvilation method, for such edge 
intiages where the numerical stability allows for it, to change the order of differentiation and 
filtrering and save 50% of the calculations for the Fourier Transform operations. Such a 
method for calculating AGVF consists of the operations: 
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w = - F'[F(f)./F(jiL-l)] (eq. 1 1). 

u.w(x,y) = w(x+l,y)-w(x-l,y) (eq. 12) and 
v.w(x,y) = w(x,y+l)-w(x,y-l) (eq. 13), 

where the symbols u_w and v_w have been chosen to emphasize that the two components of 
the AGVF field, in this case, have been calculated from a common filtering result, w. 

Since the Laplace operator has a certain geometrical width, it is not obvious that the 
components of the AGVF field are usefid or even valid all the way to the borders of the edge 
image. As a consequence, the second and the third subset of image elements, that get useful 
values of the AGVF field components may be somewhat smaller than the first subset. 

It is not necessary to use an orthogonal coordinate system x and y, neither in the microscope 
image, in the edge image nor in the AGVF field, but it makes the calculations simpler. 

In one preferred embodiment there is also performed a normalisation of the element wise 
magnitude of the AGVF field, se the field in fig. 4. Such a normalisation aims at enhancing 
the direction of the external vector force field, which is the most important information for the 
snake algorithm. 

CALCULATION OF THE ESTIMATE OF THE SOUGHT CELL CONTOUR 

In the partial operation corresponding to 87 the sought cell contour is estimated based upon a 
seed contour from a partial operation corresponding to 85 and upon the AGVF field from a 
partial operation corresponding to 86. The estimation is preferrably calculated using a snake, 
i.e. using a parametric contour model, like described in the Master* s thesis mentioned above. 
See fig. 6. 
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CALCULATION OF THE FINAL SEGMENTATION OF THE CELL NUCLEUS 

In the partial operation corresponding to 88, an estimate of the final cell nucleus contour is 
5 calculated, if needed. This can be done with a method that works based upon the estimated 
cell contour from a partial operation corresponding to 87. Se fig. 7. One possible way is to use 
a linear transformation from the three colour components of the cell image to two transformed 
colour components together with a border in the plane of the two transformed colour 
components. The transformation and the border have, preferrably, been determined earlier by 
10 experiments with manually segmented cells where the transformation as well as the border 
have been optimized in order to optimally classify cytoplasm elements and nucleus elements 
to their correct classes. The experiments may have to be reperformed for each new type of 
staining of the white blood cells. Another possible method is that the partial operation 
adaptively changes the transformation and/or the border, in order to track variations in the 
IS staining, by maximizing the spatial homogenity of the image elements that have been directed 
to each class. 

RESULTS OF THE SEGMENTATION 

20 The segmentation result from a partial operation corresponding to 89 consists of the cell 

contour, the nucleus contour and the cell image. By applying the contours, the whole cell or 
the separated cytoplasm and nucleus can be obtained as segmented images. 

Fig. 8a shows a system with an arrangement according to the invention. That system 
25 comprises a digital nucroscope 81 , which is connected to a computer system 82. The digital 
microscope 81 delivers digital nodcroscope images in gray scale or in colour to the computer 
system 82. A certain amount of digital processing of the digital images may preceed their 
delivery to the computer system 82. The computer systemet 82 may be integrated in the 
microscope 81. 

30 

The computer system 82 comprises means 83 - 89 for performing the steps that are included 
in the segmentation method described above. 
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There is included a means 83 for choosing, from a digital image, a set of image elements, 
which probably are cell nucleus elements. The output of this means is an indication of which 
elements that have been chosen in the image, see fig. 2. 



5 There is further included a means 84 for, based upon the digital cell image and upon the 
output from means 83, calculating edge values for a set of image elements that at least 
comprises the image elements that were chosen by means 83. The output of this means is an 
edge image, see fig 3* 

10 There is further included a means 85 for, based upon the output of means 83, determining a 
seed contour. The output of this means is the seed contour, see fig. 5. 

The output from means 84 is fed to a means 86, which, for the image elements that are 
included in the edge image, calculates the values of the components in an approximate GVF 
15 field. The output firom means 86 are the values of the two components in the AGVF field. 

The outputs from means 85 and 86 are fed to a means 87 that, with a snake method, estimates 
the object contour, which also is the output of the means. 

20 Using means 88, a final estimate of the cell nucleus can be determined based upon the digital 
cell image, the output from means 87 and a colour transformation and a border, the two latter 
being fixed or adaptively changing. The output from means 88 consists of the estimated cell 
nucleus. 

25 Means 89 is capable of, based upon the digital cell image and the output from means 87 and 
88, producing partial images of cell, cytoplasm and cell nucleus. The output consists of the 
partial images, which may be further analyzed in the subsequent image analysis operations 
like feature calculation and classification using artificial neural networks. 

30 The means 83-89 that are mentioned above may in general consist of the computer system, 
equipped with a computer program for performing the steps in the method. This computer 
program can be stored separately on a digital storage medium. It is however also possible to 
implement at least some of the means shown above in hardware, for example as ASIC circuits 
or as FPGA circuits. 
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It is also possible to perform calculations in a distributed way. The segmentation method does 
not have to be performed at the location of the microscope. The whole method or parts of it 
may be performed in an Internet-server, situated at a distance from the microscope. 

Fig. 9 shows a flow chart for a method 90, according to the invention, for determining a 
sought object in a digital image. 

In a first step 91 values of the elements in an edge image are determined. 

Jn a second step 92 values of the components of an AGVF field are determined. 

In a third step 93 a seed contour is determined. 

In a fourth step 94 the sought object contour is, preferrably using a snake, determined based 
upon the values of the AGVF field and the seed contour. 

The invention is useful for other applications than segmenting of white blood cells, like for 
example for segmenting of individual cell nuclei in histopathological tissue sections. 



